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ABSTRACT 


This  report  discusses  sundry  topics  in  the  theory  and  practice  of  the 
reduction  of  photometric  data.  The  topics  included  are  monochromatic 
extinction,  wide-band  extinction  with  all  second-order  terms,  transformations 
from  the  instrumental  system  to  the  standard  system,  error  analysis,  and 
observing  procedures.  In  addition,  the  approximate  computation  of  an  artificial 
satellite's  standard  magnitude  from  extremely  wide-band  data  is  discussed. 
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I. 


INTRODUCTION 


This  note  is  concerned  with  some  aspects  of  photometric  reductions 
necessary  to  obtain  accurate,  internally  consistent,  magnitudes  and 
colors.  A large  number  of  time  dependent,  temperature  dependent,  azimuthally 
dependent,  etc.,  problems  are  not  treated  here.  These  effects  are 
usually  small  (=  0™01).  Furthermore,  problems  peculiar  to  artificial 
satellites  are  not  discussed.  The  remaining  areas  considered  represent 
those  topics  capable  of  a theoretical  discussion  and  a mathematical  formulation. 
Except  for  Sec.  VI  B this  note  represents  a review  and  summary  of  the 
relevant  astronomical  literature. 

In  each  instance  the  theoretical  discussion  is  presented  first  and 
then  the  working  equations.  Several  simple  examples  are  presented  in 
great  detail  so  that  either  analytical  techniques  ot  systematic  errors 
can  be  clearly  delineated. 
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II.  NOTATION  AND  DEFINITIONS 


A(X)  * transmission  of  the  atmosphere  * dex[-0.4Xe(X) ] 
c,  C = colors  inside  the  atmosphere 
cq,  Cq  « colors  outside  the  atmosphere 

d = constant  depending  on  units,  telescope  aperture,  etc. 
D = photometric  deflection  = d 
E(X)  = emitted  energy  distribution 
F(X)  = filter  transmission  curve 


G = dJlnE(X) /d£nX 

h = height  above  earth's  surface 
H = height  of  top  of  atmosphere 
^ = scale  height  of  the  atmosphere 
m,  M = apparent  magnitudes  inside  the  atmosphere 
= -2.5  £og(D)  + constant 

mQ,  Mq  = apparent  magnitudes  outside  the  atmosphere 
n *■  -d£ne(X)/d£n(X) 


$ 


z(z')  = true  (apparent)  zenith  distance 

c,  e(X)  » atmospheric  extinction  measured  in  magnitudes 

OO  00 

X * effective  wavelength  ■ / XR(X)dX/  f R(X)dX 

° Jo  Jo 

p(h)  - atmospheric  density 


III.  PRELIMINARIES 


Before  we  can  intelligently  deal  with  the  deflection  of  the  photometer 
we  must  understand  what  happens  to  a beam  of  light  between  the  time  it 
enters  the  earth's  atmosphere  and  the  time  our  photometer  records  a 
deflection.  To  this  end  we  first  consider  which  quantities  the  deflection 
depends  on  and  then  we  describe  the  atmospheric  extinction. 

A.  The  Fundamental  Equation  of  Photometry 

The  ultimate  observational  datum  of  photoelectric  photometry 
is  the  meter  deflection  D.  The  photometer  has  converted,  according  to 
its  response  curve  P(A),  the  light  energy  falling  on  the  photocathode 
into  an  electric  current.  The  light  reaching  the  photocathode  has,  most 
recently,  passed  through  a filter  with  transmission  F(A),  a telescope 
with  transmission  T(A),  and  through  the  atmosphere  with  transmission 
A(A)  . Since  it  originally  had  E(A)  photons  at  wavelength  A it  follows 
that 

E(A)A(A)T(A)F(A)P(A)dA  = d J E(A)A(A)R(A)dA.  (1) 

The  constant  d is  fixed  once  the  unit  of  length,  the  telescope  aperture, 
etc.,  are  determined.  As  far  as  the  instrumental  system  (photometer, 
filter,  and  telescope)  is  concerned  the  effective  source  had  a distribution 

S (A)  = 


E(A)A(A). 


Astronomical  magnitude  scales  are  determined  by  the  logarthmlc 
response  of  the  human  eye  (Pogson's  Law)  and  history.  Likewise  magnitudes 
are  related  to  the  logarithm  of  the  deflection.  The  proportionality 
constant  is  -2.5  and  the  zero  of  the  scale  is  arbitrary.  Hence, 


m * -2.5  £og(D/d)  + constant. 


(2) 


Once  the  instrumental  system  is  fixed  the  additive  constant  is  too. 

B.  Air  Mass  and  Extinction 

Consider  the  passage  of  a beam  of  light  through  the  atmosphere. 
Let  the  beam  have  intensity  I (A,  h)  at  wavelength  A and  height  h above 
the  earth’s  surface.  As  the  beam  passes  through  a thickness  of  atmosphere 
ds  it  is  reduced  by  the  amount  dl  where 


dI(A,  h)  - -k<A,  h)I (A,  h)ds. 


(3) 


The  function  k(A,  h)  is  the  atmospheric  opacity.  If  the  atmosphere  is 
plane-parallel  then  ds  ■ seczdh  where  z is  the  zenith  distance  of  the 
incoming  beams  (see  Fig.  la).  With  this  expression  for  ds  we  can  solve 
for  I at  the  ground  (h  ■ 0) , viz. 


n 

I(A,  0)  - I(A,  H)exp[-secz  f k(A,  h)dh], 

A 


(4) 


M M> 





EARTH'i  surface 


where  p(h)  s 0 if  h > H,  If  R(A)  = r6(A  - A)  then  the  relation  between 
the  inside  the  atmosphere  (m)  and  outside  the  atmosphere  (m^)  monochromatic 
magnitudes  is 


m - m = -2. 5(&oge)secz  f k(A,  h)dh, 
° - h Jo 


= -1.086secz 


f k(A, 

J n 


h)dh. 


= ~e(A)secz. 


The  quantity  e(A)  is  the  monochromatic  extinction  at  wavelength  A = A.  The 
more  familiar  version  of  Eq.  (5)  is 


m = m + eX 
o 


where  X is  the  air  mass. 


The  expression  "air  mass"  for  X refers  to  the  fact  that  X is 
the  total  path  length,  measured  in  units  relative  to  the  zenith  (where  X 

f H 

= 1),  that  the  light  beam  must  pass  through.  In  general,  X “ / k(A,  h)seczdh. 
The  opacity  can  be  decomposed  into  a product  of  the  opacity  per  gram,  a, 
and  the  density,  p.  Thus, 


X « a(A)  / p(h)seczdh. 


From  Fig.  lb  we  have,  by  the  law  of  sines. 
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sinz/R  = sinz'/(R  + h), 
e e 


X « a(X)  f p(h)(l  + h/R  )dh/[cos2z’  + 2h/R  + (h/R  )2]1/2.  (9) 
J e e e 


As  long  as  the  scale  height  of  the  atmosphere,  is  small  compared  to 
Rg,  then  the  integral  in  Eq.  (9)  can  be  approximated  by 


X a a(A)p(0)secz'  [1-  (^/R  )sec2z']» 


whence. 


X = secz'[l-  (^/R  )sec2z']- 


Equation  (11)  demonstrates  two  things.  First,  it  is  the 
apparent  zenith  distance,  not  the  true  zenith  distance,  that  determines 
X.  Second,  X is  not  a linear  function  of  secz*.  However,  Eq.  (11)  is 
not  yet  accurate  because  the  curvature  of  the  ray  path  due  to  refraction 
has  been  neglected  and  Eq.  (7)  assumes  the  atmosphere  to  be  well  mixed. 
The  latter  assumption  is  demonstrably  false  for  aerosols,  ozone,  water 
vapor,  etc.,  which  contribute  to  the  opacity  by  means  other  than  Rayleigh 
scattering.  Nonetheless,  tables  of  X(z')  (cf.  Hardie  1962)  are  used 
because  there  is  no  practical  alternative.  One  should  not  use  them 
beyond  z'  - 75°  (i.e.,  X - 4). 
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IV.  MONOCHROMATIC  EXTINCTION 

The  term  monochromatic  is  used  here  to  mean  F(X)  * F6(X  - A)  where 
6(u)  is  the  Dirac  delta  function.  Although  this  is  an  idealized  case, 
it  allows  us  to  formulate  and  solve  some  problems  which  would  be  too 
complex  to  handle  otherwise.  The  treatment  of  monochromatic  extinction 

also  lays  the  groundwork  for  the  more  complex  wide-band  case.  Hence,  it 

. 

introduces  the  methods  (but  not  reduction  procedures)  necessary  in  the 
ab  initio  construction  of  photometric  systems.  The  best  photometric 
system  for  artificial  satellites  will  be  discussed  elsewhere. 

A.  The  High-Low  Method 

This  represents  the  simplest  procedure  and  does  not  assume 
that  the  extra-atmospheric  magnitude  is  known.  Moreover,  it  is  necessary 
to  observe  only  a single  star.  Let  the  star  be  observed  at  air  masses  X 
■ X^,  X£  with  corresponding  magnitudes  m =■  m^,  m^.  Since, 

m^  = mQ  + eXj,  j =*  1,  2,  (12) 

clearly 


e » (n^  - m2)/(X1  - X2). 


(13) 


If  | X^  — X2 | is  small  then  |m^  - m2|  will  be 
with  which  e is  known  would  be  poor.  On  the 
in  the  air  masses  implies  several  hours  have 


also.  Hence,  the  accuracy 
other  hand,  a large  difference 
elapsed  between  the  two  observations. 
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Therefore,  problems  associated  with  Instrumental  drift,  time  varying 
extinction,  etc.,  are  encountered.  Ignoring  these  difficulties,  it  Is 
still  not  wise  to  maximize  [x^  - X2(  because  we  want  to  minimize  0£,  the 
estimate  for  the  standard  deviation  of  the  extinction. 

Let  us  expand  upon  this  point  further.  It  is  a common  (and 
incorrect)  practice  to  assign  equal  weights  to  all  observations  when 
determining  e.  However,  it  is  known  that  the  standard  error  of  a measurement 
depends  on  the  air  mass.  Stock  (1968)  has  found  a quadratic  air  mass 
dependence.  Thus,  if  0 ^ is  the  standard  error  of  a measurement  made  at 
the  zenith  (X  = 1) , the  weight  for  a measurement  made  at  any  other  air 
mass  will  be 


9 2—2 

w « 1/a  = (a.X  ) . 


Consider  the  case  X^  = 1,  * X.  Then, 


°e  " °l(x4  + 1)/(X  " 1)2, 


The  extinction  will  be  determined  most  accurately  when  0£  is  a minimum. 

This  implies  X^  - 2X^  -1  * 0 or  X = 2.10692.  Thus,  if  X * secz,  the  off 

zenith  position  corresponds  to  z * 61°40'  and  ■ 4.111a^.  The  effect 

2 

of  choosing  another  air  mass  can  be  seen  in  Fig.  2 where  *-s 

plotted  versus  X. 
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As  large  air  mass  observations  carry  low  weight.  It  would  be 
natural  to  spend  more  time  observing  there.  In  order  to  analyze  this, 
and  optimize  telescope  time,  suppose  N observations  are  available  with 
the  fraction  f to  be  made  at  the  X > 1 position.  In  this  case 

a l * 02[X4/f  + 1/(1  - f)]/[N(X  - 1)21,  (16) 

and  we  want  to  minimize  the  estimate  for  the  standard  deviation  of  the 
extinction  with  respect  to  X and  f.  The  result  is  X * [(2  +^2)/ 

(2  -^2))1/2  * 2.41421  (z  = 65°32' ) and  f = (2  +^2)/4  - 85.4  percent . 

The  value  of  0£  = 4. 828a ^ See  Young  (1974)  for  a fuller  discussion 

of  this  case. 

These  two  examples  suffice  to  demonstrate  the  necessity  of 
ascertaining  the  air  mass  dependent  errors  when  computing  the  extinction 
in  even  the  simplest  cases.  Furthermore,  an  indication  of  the  division 
of  telescope  time  can  be  obtained. 

The  next  step  upward  in  complexity  is  to  use  a pair  of  stars 
and  observe  each  of  them  twice.  Let  the  stars  have  magnitudes  m,  M and 
be  observed  at  air  masses  x,  X.  The  first  set  of  observations  occurs 
when  the  lower  case  star  is  near  the  zenith  (m  = m^,  x = x^)  and  the 
upper  case  star  is  several  hours  east  of  the  meridian  (M  = M^,  X - X^) . 
The  second  set  of  observations  takes  place  when  the  upper  case  star  is 
near  the  zenith  (M2,  X2)  and,  consequently,  the  lower  case  star  is 
several  hours  west  of  the  meridian  (m2,  X2).  Then, 


* 


f.1: 

i 


l 


E ■ [(n^  - - (n>2  - M2)/[(x1  - Xx)  - (x2  - X2).  (17) 

This  procedure  has  the  advantage  of  eliminating  time  dependent  effects 
in  the  telescope-photometer  system  because  the  individual  observations 
within  a set  are  nearly  simultaneous.  Also,  the  extra-atmospheric 
magnitudes  mQ,  MQ  need  not  be  known. 

Once  more  than  two  stars  are  used  to  determine  £ an  algebraic 
solution  is  no  longer  possible  and  least  squares  must  be  used. 

B.  Simultaneous  Standard  Stars 

Let  us  suppose  that  several  stars  are  chosen  to  be  standards 
and  observed  over  enough  nights  to  determine  good  values  for  m^.  Good 
means  that  their  extra-atmospheric  magnitudes  have  been  calculated 
sufficiently  accurately  to  rule  out  the  possibility  of  variability. 

Then,  to  determine  the  extinction,  we  observe  several  of  these  stars 
over  a range  of  air  masses.  We  now  want  to  minimize 
£w[m  - mQ  - eX]2/£w 

with  respect  to  £.  The  w's  are  the  air  mass  dependent  weights. 

However,  if  the  extra-atmospheric  magnitudes  are  unknown, 
iterative  procedures  are  available  (e.g..  Weaver  1952)  to  simultaneously 
determine  both  £ and  {mo} • For  this  purpose  observations  are  performed 
over  several  nights  on  the  same  set  of  stars.  Of  course,  unless  the 
absolute  sensitivity  of  the  equipment  can  be  established  by  independent 
means,  nightly  zero  point  corrections  must  be  included.  Thus,  for  the 
j'th  observation  on  the  k'th  night  of  a star 
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(18) 


IT 


F 

r ' 

|e  ; 


4 


"jk 


m + e.  X,.  + z,  . 
o k jk  k 


C.  Monochromatic  Colors 

In  astronomical  photometry  a color  (e.g.,  color  index)  is 
the  difference  between  two  magnitudes.  That  is,  if  the  two  monochromatic 
magnitudes  of  a star  are  m,  M then  c * m - M is  a color.  From  Eq.  (4) 
and  the  definition  of  c it  follows  that 


1 


. 


3 

- 


I 

I 


c = c - ex. 
o c 


(19) 


Hence,  a high- low  method  can  be  used  to  determine  e^.  We  note  that 
sensitivity  drift  does  not  affect  Usually  m refers  to  a shorter 
wavelength  than  M. 
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V.  WIDE-BAND  EXTINCTION 

When  the  instrumental  response  function,  R(X),has  a bandwidth  (> 

o 

100  A)  sufficiently  large  that  the  extinction  varies  across  it,  the 

determination  of  the  effective  extinction  is  more  complicated.  It  is 

obvious  from  the  appearance  of  the  sun  at  noon  and  sunset  that  e(X  . ) 

blue 

51  e^recj)*  Hence,  the  use  of  a constant  E in  Eq.  (6)  introduces  two 
different  errors.  The  first  is, because  de/dX  # 0,a  blue  star  suffers 
more  extinction  than  a red  one  (Stebbins,  et  al.  1950,  Eggen  1950).  The 
second  is, that  since  X ^ secz',  as  a star's  altitude  decreases  it  becomes 
redder  and  the  effective  extinction  becomes  smaller  (Forbes  1842).  To 
ascertain  the  importance  of  these  complications  we  follow  King  (1952) 
and  correctly  calculate,  through  terms  of  the  second-order,  the  analog 
of  Eq.  (6) . 

A.  Theory 

The  photometer  deflection,  D,  is  given  by 

00 

D - d f S(X)R(X)dX.  (20) 

o 


Since  R(X)  peaks  at  some  wavelength  near  X = X let  us  expand  S about 

o 

X . 

o 


D * d f [ S(X  ) + (X  - X )S ' (X  ) + (X  - 
J . o o o 


V s"(xo>/2]R(x)dx» 


d[S  + WX*S 
o o o 


dE  A [1  + WX  S 
o o o o 


onV  f 


R(X)dX, 


R(X)dX, 


(21) 
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where  SQ  = S(Ao)  etc.  The  under  the  atmospheric  magnitude  m = -2.5  £og(D/d) 
plus  constant  and  the  magnitude  outside  the  atmosphere,  mQ,  is 
given  by  a similar  expression  with  A(A)  * 1.  Hence, 


m - m - 2.5  AogA  - 2.5  Jlog{[l  + WA2S"/2S  ]/ [1  + WA2e"/2E  ]}, 


o o o o o 

= m + e X - 1.25WA2(Jloge)[s"/S  - e’  /E  ], 
o o o ° 1 o o oo 


o o o 


where  e X * -2.5  £.ogA  . Since  S = EA  this  can  be  rewritten  as 
o o 


m - m - e X[1  - n GW  + n (n  + l)W/2  - WA  n'/2] 
oo  o o o oo 

- WHnlO[n  e X]2/5. 
o o 


The  various  terms  on  the  right  hand  side  of  Eq.  (23)  represent 

(i)  the  monochromatic  extinction  evaluated  at  the  effective  wavelength, 

(ii)  a correction  which  depends  on  the  color  of  the  star  through  the 
gradient  of  its  energy  distribution,  (iii)  a second  correction,  present 
for  E(A)  * constant,  which  depends  on  the  color  of  the  atmosphere  through 
the  gradient  of  its  absorption,  (lv)  a third  air  mass  dependent  correction 
which  depends  on  the  square  of  the  atmospheric  color,  and  (v)  a non- 
linear term  in  eX  (e.g.,  the  Forbes  effect).  In  order  to  establish 
(heuristlcally)  that  G depends  on  the  stellar  color  we  approximate  E(A) 


by  the  black  body  distribution. 

The  black  body  spectrum  is 


E(A)  = CjA  /[exp(c2/AT)  - 1], 
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where  and  are  the  first  and  second  radiation  constants  and  T is 
the  absolute  temperature.  A straightforward  calculation  yields  (y  = 

c2/XoT)’ 


G - -5  + yey/(ey  - 1). 


For  the  majority  of  the  stars  we  observe  the  effective  wavelength  will 
be  near  the  wavelength  satisfying  the  Wien  displacement  law.  Thus,  y - 
5,  ey  >>1,  and 


G “ -5  + y. 


The  star’s  color,  c , is  given  by  c = m (A ’ ) - m (A").  Thus, 
o o o o 


-1.25W(£oge)[(A,)2E"(A’)/E(A’)  - (A")2E"(A")/E(A") ] 
-1.25WOLoge)[-12(y’  - y")  + (y’)2  - (y")2] 

- 30W (loge) yQAA / A q , 


where  A’  = A + AA,  A"  ■ A - AA . Thus,  comparing  Eqs.  (26,  27) 
o o 


G * a + be  , 
o 


since  y , AA,  and  W are  fixed, 
■'o’ 
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_2  in 

For  n ■ 4 (Rayleigh  scattering),  W - 10  and  e - 0.25.  The 


Forbes  effect  term  varies  by  only  O^Ol  as  X varies  from  1 to  3.  If  e * 
0?5  then  over  the  same  range  of  air  mass  it  changes  by  0?15.  Hence,  on 
a good  night  it  represents  no  problem  but  on  a poor  night  it  results  in 
a systematic  error  in  mQ  (or  alternatively  eQ) . The  n^  term  is  difficult 
to  evaluate  since  it  depends  on  those  components  of  the  atmosphere  whose 
composition  is  the  most  variable.  Thus,  its  systematic  neglect  will 
reappear  in  the  large  standard  deviation  of  for  stars  of  widely 
different  color. 

B.  Practice 

It  is  standard  astronomical  practice  to  replace  the  exact 
higher  order  terms  in  Eq.  (23)  with  observable  quantities.  Thus,  one 
replaces  Eqs.  (6,  19)  by 


m - m + (e'  + e"c  )X,  (29a) 

o o 

c - c + (e*  + e"c  )X,  (29b) 

o c co 

where  e',  are  referred  to  as  the  principal  (or  first-order)  extinction 

coefficients  and  e",  e"  are  referred  to  as  the  second-order  extinction 

c 

coefficients.  Sometimes  the  color  index  under  the  atmosphere,  c,  is 
used  in  place  of  cq  in  Eqs.  (29).  From  Eq.  (23)  we  see  that  both  the 
stellar  energy  distribution  and  the  atmospheric  absorption  affect  the 
effective  extinction  coefficients.  Hence,  if  cq  Is  used  in  Eqs.  (29) 
only  the  stellar  component  is  represented  whereas  if  c is  used  both  the 
stellar  and  atmospheric  components  are  included. 
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C.  Determining  the  Extinction 


Consider  two  stars,  of  different  color,  close  together  on  the 
sky.  Then  if  one  star  has  lower  case  magnitudes  and  colors  while  the 
other  star  has  upper  case  ones,  we  can  write 


m -M  - m - M - e"(c  - C)X, 
o o 

(30a) 

c -C  ■ c - C - e"(c  - C)X. 
o o c 

(30b) 

Several  measurements  performed  at  different  air  masses  allow  the  least 

squares  determination  of  the  second-order  terms  (with  m , C , etc., 

o o 

unknown).  Once  these  have  been  determined,  an  optical  double  suffices  to 
compute  the  first-order  terms  by  the  high- low  method.  Finally,  since 
the  second-order  coefficients  vary  (night  to  night)  less  than  the  principal 
coefficients  do,  an  analog  of  the  high-low  method  will  quickly  yield  the 
first-order  coefficients  once  the  second-order  ones  are  fixed. 
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VI.  MAGNITUDE  AND  COLOR  TRANSFORMATIONS 


A.  Theory  and  Practice 

When  one  attempts  to  do  photometry  on  a standard  system  one  Is 
faced  with  the  problem  of  reducing  one's  own  magnitudes  and  colors  to 
the  standard  system.  This  is  best  accomplished  by  matching  R(A) , reduction 
procedures,  observing  procedures,  etc.,  as  closely  as  possible.  In  the 
end,  the  observer  has  a magnitude  m and  color  c in  his  system  that 
should  be  close  to  the  corresponding  standard  quantities,  M and  C.  From 
the  analysis  in  Sec.  V it  can  easily  be  shown  that  (for  black  bodies) 
the  following  linear  transformations  suffice: 


M = m + ac  + $ = m + AC  + B, 
C = yc  + 6. 


(31a) 

(31b) 


Clearly  6(B)  and  6 are  zero-point  constants  and  a(A)  and  y are  scale 
factors.  The  values  of  y and  a relative  to  unity  and  zero  reflect  the 
sizes  of  the  two  bandwidths  and  the  relative  position  of  their  effective 
wavelengths.  Thus,  if  y < 1 the  standard  system  has  a narrower  bandwidth 
than  the  instrumental  system  while  if  a < 0 the  effective  wavelength  of 
the  standard  system  is  redder  than  that  of  the  instrumental  system. 

A mismatch  in  R(X)  that  is  easy  to  deal  with  is  when  W remains 


constant  but  Aq  changes.  Let  R(A)  be  the  instrumental  response  function 
of  the  standard  system  with  effective  wavelength  Aq.  Then  if  r(A)  is 
the  observer's  response  function  with  effective  wavelength  Aq, 


D = d Jq  S(A)r(A)dA 

= d / [S  (A  ) + (A  - A )S ' (A  ) ]r  (A)dA 
Jo  ° o o 

r°° 

= dS  [1  + (A  - A )S’/2S  ] / r (A)dA 

O O O O O Jo  oo 

= dS  {1  + (A  - A ) [G  + 0.4n  e X«hlO]}  / r(X)dA. 
o o o o o Jo 


Hence, 


M=m-2.5(A  - A )[G  + 0.4n  e XinlO], 

o o o o 


and  we  already  know  [cf.  Eq.  (28)]  G Is  a function  of  the  color.  The 

additional,  air  mass  dependent,  term  varies  over  Q^ll  as  X varies  from  1 
oo 

to  3 for  a 100  A mismatch  at  5000  A with  Rayleigh  scattering  and  0.5 
extinction  per  air  mass.  This  can  dominate  the  color  term  and  furthermore, 
it  is  clearly  systematic  in  X. 

Finally,  we  consider  the  correct  least  squares  determination 
of  the  constants  in  Eqs.  (31).  We  assume  that  the  standard  values  are 
without  error.  Then  one  should  use  the  regression  line  of  c on  C to 
determine  y and  6.  Hence, 


2 2 
Y = ac/cov(c,  C)  j c0v(c,  C )/ac. 


(34a) 


where  0^(0^)  is  the  variance  of  C(c)  in  the  sample  and  cov(c,  C)  is  the 
covariance  of  c and  C.  The  value  of  6 is  calculated  from 
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5 - <C>  - y<c> 


(34b) 


where  angular  brackets  denote  a sample  mean.  When  computing  A and  B It 
must  be  recognized  that,  since  C depends  M,  the  errors  In  M and  C are 
correlated.  Hence, 


[0^  - cov(m,  M) ]/cov(M,  C)  + cov(M  - m,  O/a^, 


(35a) 


while 


B - <M>  - <m>  - A<C> . (35b) 

It  would  be  simple  to  obtain  estimates  for  the  standard  deviations 
of  these  constants  and,  thence,  formal  estimates  for  the  errors  in  M and 
C.  It  Is,  however,  more  meaningful  to  observe  the  same  stars  on  several 
nights,  reduce  each  night  separately,  and  then  calculate  the  internal 
variance  of  M and  C.  These  quantities  are  frequently  used  as  estimates 
of  the  external  errors  of  the  transformations. 

B.  S-20  Transformations 

The  principal  artificial  satellite  detector  for  GEODSS  is  an 
S-20  photocathode.  To  maximize  the  signal-to-noise  ratio,  and  hence  the 
ability  to  confidently  detect  faint  satellites,  one  does  not  want  to  use 
filters  in  the  light  path.  Nonetheless,  for  some  considerations,  it 
would  be  advantageous  if  at  least  a magnitude  could  be  reported,  in  a 
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*sw* 


; 


TABLE  1 


■ 


STARS  WITH  NEARLY  SOLAR  TYPE  SPECTRAL  DISTRIBUTIONS 


B.  S.  # 

100W 

V 

100W_ 

D 

B.  S.  # 

100W 

V 

100W] 

77 

30 

33 

4707 

24 

68 

98 

13 

10 

4708 

34 

39 

219 

32 

29 

4785 

26 

20 

244 

32 

44 

4845 

40 

42 

321 

38 

50 

4833 

23 

25 

370 

24 

23 

4979 

51 

50 

424 

29 

44 

4983 

24 

31 

458 

41 

52 

5019 

44 

50 

483 

14 

9 

5235 

23 

44 

618 

57 

45 

5260 

22 

31 

641 

39 

38 

5304 

44 

55 

788 

14 

22 

5338 

44 

58 

937 

14 

15 

5384 

20 

14 

996 

29 

32 

5409 

46 

55 

1017 

36 

87 

5868 

16 

17 

1101 

24 

32 

5914 

34 

42 

1242 

46 

73 

5986 

47 

60 

1674 

50 

63 

6098 

53 

37 

1729 

11 

6 

6978 

35 

39 

2047 

19 

23 

7172 

49 

61 

2693 

49 

49 

7503 

17 

20 

2906 

38 

52 

7504 

21 

18 

3064 

24 

19 

7796 

50 

48 

3391 

16 

13 

7955 

36 

47 

3591 

43 

43 

8131 

27 

66 

3871 

30 

75 

8170 

42 

54 

3881 

12 

13 

8334 

41 

28 

4112 

49 

62 

8729 

20 

22 

4540 

31 

41 

8737 

12 

11 

8817 

36 

35 

8905 

14 

8 

24 

5 


fitting,  a weighting  « 1/ (monotonic  function  of  W ) would  be  appropriate. 

D 


A standard  G2  V star  has  W = O^IS. 

D 


There  is  no  doubt  about  the  success  of  this  procedure.  However, 

■>-2 


unless  max  ( |a  | , |av | , | ar | , |a^ | ) = 10  it  will  not  yield  useful 


information.  Finally,  if  V were  desired  one  would  construct  a analogously. 


Table  1 also  includes  stars  for  which  <_  0.5. 


Since  Eq.  (36)  is  an  extension  of  a relationship  for  black 
bodies,  a more  fundamental  approach  may  yield  better  results.  In  terms 
of  extra-atmospheric  values, 
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L 


E(X)P20(X)T(A)dX, 


f E(X)Fb(X)P2Q(X)T(X)dX,  etc. 

Jo 


(38) 

(39) 


Thus,  we  can  regard  D2q  as  given  by  a five-point  quadrature  formula, 


D0  = A D + A,  D,  + A D + A D + A.D.  . 
20  uu  bb  vv  rr  ii 


(AO) 


This  implies 


B = m2g  + bz  - 2.5  J!.og{budex[-0.4(U  - B)  ] + b^ 


+ b dex[-0.4(V  - B) ] + b dex[-0.4(R  - B) ] 


+ b1dex[-0.4(I  - B)]}. 


(41) 


The  procedure  outlined  above  would  again  be  followed  but  the  least 
square  problem  is  more  complicated.  Moreover,  the  systematic  accuracy 
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